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Abstract 

The direct detection of gravitational waves will provide valuable astrophysical information about 
many celestial objects. Also, it will be an important test to general relativity and other theories 
of gravitation. The gravitational wave detector SCHENBERG has recently undergone its first test 
run. It is expected to have its first scientific run soon. In this work the data analysis system of this 
spherical, resonant mass detector is tested through the simulation of the detection of gravitational 
waves generated during the inspiralling phase of a binary system. It is shown from the simulated 
data that it is not necessary to have all six transducers operational in order to determine the 
source's direction and the wave's amplitudes. 
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I. INTRODUCTION 



As predicted by the theory of relativity and other theories of gravitation, time-dependent 
gravtafoual forces are expected to propagate m spacetime in the form of waves Q.^Strong 
gravitational waves (g.w.) are expected to be generated by astrophysical objects [2']. For 
instance, two very massive stars orbiting each other would emit such waves. In particular, 
when they are coalescing they emit waves in a large frequency range. Another example is 
given by a black hole that rings down, also emitting in different wavelengths jslj. 

There is the belief within the international community that works with g.w. detection 
that the first direct detection of gravitation radiation from an astrophysical source will 
become a reality in the near future. This will open a new window to the universe, bringing 
new information about known objects and about fairly unknown things, like dark matter. 
In order to extract such information from the huge amount of data that is expected to be 
generated from the detectors' outputs (which has already started) a lot of work will be 
demanded in the field of data analysis. 

From the analysis of the detected waves emitted by such sources important information is 
expected. The very first direct detection will provide a test for one of the predictions of the 



theory of general relativity. Then, continuous observation will al 



ow testing other theories 



of gravitation Q,y,B, besides initiating gravitational astronomy flfl. 

In order to make astrophysical observations the following parameters that characterize 
the g.w. are needed: the amplitudes of the two states of polarization of the wave as functions 
of time and hy,{t)) and the source's direction in the sky (given by the angles G and 

$) 9|. Therefore, gravitational wave observatories must be able to detect four independent 
observable in order to allow for gravitational wave astronomy to start. 

The basis for the existence of resonant-mass gravitational wave detectors like 



SCHENBERG 



10| is the fact that solid bodies are distorted in the presence of g.w. due 



to the changes in spacetime. In principle, the spherical geometry of SCHENBERG's an- 
tenna implies no preferred direction of observation (i.e., it is omnidirectional) and the all 
observable needed for gravitational astronomy can be obtained from only one spherical de- 
tector appropriately equipped 0, [ill- When fully operational, this detector will be able not 
only to acknowledge the presence of a g.w. within its bandwidth: it will be able to inform 
the direction of a source in sky, the wave's amplitude and the polarization components in 



2 



the detector bandwidth - one only antenna working as a gravitational wave observatory in a 
bandwidth between 3000Hz and 3400Hz, sensitive to displacements around 10~^°m. To this 
end six transducers will continuously monitor radial displacements of the antenna's surface. 
SCHENBERG is installed at the Physics Institute of the University of Sao Paulo (Sao Paulo 
city, Brazil) and has undergone its first test run in September 8, 2006, with three transducers 
tuned to the above band. 

There are mathematical models for this detector for the case that all these six transducers 
are operational. Such models have been investigated two situations: one in which the 
transducers are perfectly uncouj 3led [nl, Q, Q and another in which the transducers are 
somehow coupled to each other |14l]. For both cases it is possible to solve the inverse problem 
and obtain the four astrophysical parameters needed if one assume General Relativity as 
the correct theory of gravitation. 

This work presents the results of an investigation 15| carried on by the data analysis 
group within the GRAVITON project IGj (the one SCHENBERG is part of) aimed at 
developing a model of the detector with less than six transducers. This is an important 
issue, since occasionally not all six resonators may be operational simultaneously. In this 
case there is a break in the convenient symmetry among the transducers. Two approaches 
were under investigation: one considered the model already developed for 6 transducers [l7| 
and simply reduced their number; the other, described here, considered the fewer transducers 
as independent devices and redesigned the mathematical model. It was found that the 
system sphere plus transducers is equivalent to a system of independent bar detectors with 
center-of-masses located at the sphere's center-of-mass ("cyclical" bars) [7]. It was possible 
to retrieve the relevant astrophysical information about a coalescing binary neutron star 
system both from the sphere coupled to six transducers and from a system considering the 
sphere's outputs as the outputs of transducers coupled to cyclical bar detectors (which will 
be called the "method of independent bars"). These results are detailed as follows. 



II. METHOD OF INDEPENDENT BARS 

n 

It has already been shown theoretically by Magalhaes et al.[7] that any array of n cyclical 
bars will respond to a gravitational wave in a way similar as a non-noisy, high Q spherical 
antenna monitored by n transducers. In both cases, the projection of a g.w. onto a 
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detector can be expressed in terms of the internal product between two symmetric, traceless 
tensors Q]: 

3 

R' = J2 ^^^ A- ■ (1) 

In this expression are the outputs of the transducers (L=l,2,...,n, where n is the number 
of transducers positioned on the sphere's surface or, equivalently, the number of cychcal 
bar detectors under consideration). The tensor W contains only information concerning the 
gravitational wave, while carries only information concerning the detector. For a bar 
detector, if 6l is the angle the detection point (where the L-th transducer is located) makes 
with the local zenith and 0l is the azimuthal angle then one can show that 



i i sin^ sin 2(^2, ^ sin 20^;, cos 0l 
i sin^ 6*^ sin 2(/)2, sin^ 0^ sin^ 0/, — i i sin26'L sin0L 



(2) 



2 Olll W Ij Olll ^^Jj Olll L/ 0±11 ^Jj 3 2 

i sin 26*^ COS ^ sin 20^;, sin cos^Oi 

It is the knowledge of the wave's tensor, W ^ that allows one to determine all astrophysical 
parameters, (/i+(t), h-^{t)) and (6, $), as was shown in [3]. In order to determine Wij 
one notices that, from general relativity, this tensor has six independent terms. It is then 
necessary to have at least six equations to determine these six unknowns. One equation is 
the traceless condition 



Wii + W22 + 1^33 = 0, 



(3) 



the other five can be given by 



(4) 



^2^2^22 + 2^2^3/123 + 2L'|'3^33 = 2i?^. 

with L = 1, ...,5. As usual, at least five cyclical bars with one transducer each are needed. 
The six equations given by ([3]) and (jlj) can be written in matrix form as 

^33 



Dl 2Dl, 2Dl, Dl^ 2Dl, Dl 



Df, 2Df, 2Df, Dl, 2Dl, D: 
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FIG. 1: Transducers' outputs from the simulation using as input a g.w. emitted by a coalescing 
NS-NS system. 

or DW = R. As long as the matrix D is square and non-singular, from the experimental 
knowledge of and one can obtain the wave's tensor W through the relation 



W = D-^R. (6) 

In order to make sure that the use of a tensor for bars, ([2]), would yield the solution 
for the inverse problem using R^ derived from the data of a spherical detector, a test was 
performed. The transducers were positioned as planned in SCHENBERG, located according 
to the pentagonal faces of a truncated icosahedron: 

[(e, (f))] = [(79.1877°, 120°), (79.1877°, 240°), (79.1877°, 0°), 

(7) 

(37.3774°, 300°), (37.3774°, 180°), (37.3774°, 60°)]. 

For future use in the calculation of ([6]), in the construction of the matrix D only the trans- 
ducers located in the five first positions were used. 

The values of R^ were obtained using the software "Schenberg Simulator 2.0.0" [l3], ini- 
tially in the absence of noises. The mathematical model of the detector, which is the basis 
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for the software, can be found, for instance, in 



In the software are used six two-mode 



transducers located on the antenna's surface, distributed according to the angles listed 
above, operating simultaneously. This software has proven effective in solving the inverse 
problemflT]. The outputs of five of these transducers (the same used to construct the matrix 
D) will be used as the components of . 

If the equivalence between the spherical detectors and the independent (cyclical) bars is 
correct, then the solution of the inverse problem with the bars should yield the astrophysical 
parameters of a source introduced as the input of the software that simulates SCHENBERG's 
response. As will be shown, this equivalence does exist. 

The simulation that is going to be analyzed here was run with the source represented by 
a waveform resultant from a coalescing neutron star (NS) binary system [3], incident on 
the detector according to the angles O = 60° and $ = 45° relative to the same system that 
defines the transducers' positions. With the input from this source the software yielded the 
outputs illustrated in Figure [H 

Using the first all the six outputs in eq. ([6]) the values of W were found and, 
consequently 71] , the source's direction was retrieved. In this particular case considered we 
have a probability of 68.27% that the direction is 6 = 60° ± 2° and $ = 44.4° ± 1.4°. For 
the case of five outputs in eq. ([6]) the result of the source's direction is = 60° ± 2° and 
$ = 44.2° ± 1.6°. It was verified that the result was in perfect agreement with the source's 
input angles. The wave's amplitudes were also determined, as shown in Figure [2l 

Different combinations of five transducers were used, yielding the same results as expected 



ll| . It became clear that the method of independent bars using five outputs was as effective 



in solving the inverse problem as the conventional method, which considers all six transducers 
simultaneously. 

The procedure was repeated in the presence of Gaussian noises and the results from both 
approaches coincided again. See Figure [3] for the case when the signal-to-noise ratio, SNR, 
is approximately 18. In this figure two regions of the sky are marked because there is a 
natural ambiguity in the determination of the direction by only one detector; if more than 
one detector is used then that there is a time delay between detections and one can determine 
whether the wave came from up or down. The wider the region marked, the higher the SNR. 

The method of independent bars was also tested using the output of only four transducers 
to solve the inverse problem in the presence of Gaussian noise for the case when SNR is the 
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FIG. 2: Amplitudes of the g.w. emitted by the simulated coalescing NS-NS system. These plots 
were obtained from the solution of the inverse problem with the method of independent bars using 
five outputs in the absence of noise. 
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FIG. 3: Source's direction as determined by the method of independent bars using five bars in 
the presence of Gaussian noise with SNR~ 18. The different degrees of gray are related to the 
differences in the energy absorbed by the sphere. 

approximately 93. The preliminary results show that even using less transducers all the 
four astrophysical parameters were determined from the method. The result of the source's 
direction in this case is = 57° ± 5° and $ = 42° ± 6°. This result take obtained using the 
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following transducers: 

[(e, <f))] = [(79.1877°, 120°), (79.1877°, 240°), (79.1877°, 0°), 

(8) 

(37.3774°, 60°)]. 

However, in this case the computational time increases considerable. Also, the precision of 



the results vary with the choice of the transducers, confirming previous investigation [ll|. 



III. SUMMARY 



Presently, two spherical detectors exist (SCHENBERG and MiniGrail), and they are 
expected to begin scientific runs soon. This work presents results of a method that allows 
determination of astrophysical parameters from the data of a spherical g.w. detector using 
6, 5 or 4 transducers. This is an original method, since in the literature only 6 transducers 
were considered so far. 

The research will continue with the solution of the inverse problem using the method of 
independent bars with the output of 4 and 3 transducers, both assuming a noiseless and a 
noisy detector. Circular and elliptical polarizations should be considered, in addition to the 
linear ones. The detailed presentation of these results will be published in a longer paper 
that is now in preparation. 

Also, it is intended to fully investigate SCHENBERG's directivity pattern, including 
the simulation of real data from signals arriving from different directions in the sky. This 
will allow statistical estimates of the actual detector performance based on the theoretical 
principles exposed here. Such an statistical analysis is essential since the less transducers 
one has, the more important, if not crucial, is to assess the significance of the detection. 
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